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Abstract 

We demonstrate how HQET and the Step Scaling Method for B-physics, pi- 
oneered by the Tor Vergata group, can be combined to reach a further improved 
precision. The observables considered are the mass of the b-quark and the B s - 
meson decay constant. The demonstration is carried out in quenched lattice QCD. 
We start from a small volume, where one can use a standard 0(a)-improved rel- 
ativistic action for the b-quark, and compute two step scaling functions which 
relate the observables to the large volume ones. In all steps we extrapolate to 
the continuum limit, separately in HQET and in QCD for masses below m^. The 
physical point is then reached by an interpolation of the continuum results 
in 1/m. The essential, expected and verified, feature is that the step scaling 
fuctions have a weak mass-dependence resulting in an easy interpolation to the 
physical point. With r$ — 0.5 fm and the experimental B s and K masses as in- 
put, we find F& s — 191(6) MeV and the renormalization group invariant mass 
Mb = 6.88(10) GeV, translating into m h (m h ) = 4.42(6) GeV in the MS scheme. 
This approach seems very promising for full QCD. 



1 Introduction 



It has long been realized that B-meson decays and mixing have a significant potential 
for the search for physics beyond the Standard Model of particle physics. Unfortunately, 
the comparison of experimental results from BaBar and Belle to the Standard Model 
has not yet revealed such effects. An even higher precision in both future experiments 
and the corresponding "predictions" of the theory is required if we want to get hints for 
new particles or interactions in this wayQ 

The most promising method for the computation of QCD matrix elements with at 
most one hadron in initial and final states is lattice QCD. Contrary to what is some- 
times reported, it is, however, a very non-trivial task to achieve precisions at the (few) 
percent level, keeping all systematic uncertainties under control. This is particularly 
so in B-physics, where the difficulty of simulating light quarks with masses that make 
contact to the regime where chiral perturbation theory is applicable meets the addi- 
tional requirement of correctly describing the physics of the heavy b-quark. The former 
requires lattices of a large enough physical size, say 2 — 3fm across and the latter a 
small lattice spacing, a, or the control of an effective theory (see |1|2|3|4|5] for more 
detailed accounts of the difficulties and recent progress). In this letter we exclusively 
discuss a method to cope with the discretization errors associated with the heavy quark 
dynamics. The light quark is simply taken to be the strange quark and for the purpose 
of testing the methodology we work in the quenched approximation. The light (dy- 
namical) quark simulations are an entirely separate issue, where fortunately significant 
progress has recently been made |6)7)8)9|10|ll|12|13)14j . 

The basic idea of the approach investigated here is the fact that b-quarks can 
be simulated (quite) straightforwardly in a space-time volume with a linear extent of 
Lq = 0(0.5 fm) |15|16j . In such a volume the lattice spacing can be chosen small enough 
such that observables can be computed with a relativistic action for the heavy quark. 
The continuum limit is reachable by a short, controlled, extrapolation. Starting from 
this simple idea, two different roads have been taken in the past |17|18|19] and a third 
one has recently been explored [20] . 

In the first the (continuum) observables in the small volume serve to determine the 
parameters of HQET non-perturbatively and then the physical (large volume) matrix 
elements are computed in this effective theory. By including 1/mb-corrections a good 
overall precision is attainable |19j . 

In the second way, one remains in the relativistic theory, and computes the finite 
size effects of the observables iteratively (Lq — > L\ = sLq — > L2 = s 2 Lq . . .). As one 

1 In fact the situation on the theory side is not sufficiently clear to exclude that experiments have 
found new physics already. The point is that hadronic matrix elements of B-mesons are difficult to 
compute. It could thus be that some matrix element (B-factor or other) which has been extracted from 
fits to the unitarity triangle is actually in disagreement with the true matrix elements in QCD. Improved 
determinations of these matrix elements are hence of interest even without an increase of precision of 
the experiments. 
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increases the volume also the lattice spacing is increased and one has to reduce the 
mass, mji, of the actually simulated quark to remain with am^ <C 1. The physical mass 
of the b-quark is then reached by an extrapolation. 

Here we demonstrate how the two approaches can be combined by constraining 
the extrapolation to the physical quark mass with calculations in the effective theory; 
extrapolations are turned into interpolations and an even higher precision as well as 
confidence is reached. 



2 Strategy 

We are interested in computing an observable O, which, in addition to the light quark 
masses, depends on the mass, m^, of a heavy quark. Its exact definition will be men- 
tioned when it becomes relevant. In a Monte Carlo computation the observable depends 
in addition on the linear extent L of the simulated space-time volume. This finite size 
effect is negligible when L is large enough, which we here assume to be the case for 
L > Ljy. Following |21|22j . we express O as a product of factors, 

0(m h , Loo) = 0(m h , L ) — — • • • ■ (2.1) 

Here Lq is chosen small enough such that with an affordable effort lattices with a spac- 
ing o< 1/ m-h can be used and the continuum limit can be reached by an extrapolation 
of O computed with a relativistic 0(a)-improved action. For the b-quark this means 
that a ~ 0.012 fm can be used. In a small volume the details of the topology, boundary 
conditions and the exact choice of observables are relevant. We here note only that 
choosing Schrodinger functional boundary conditions makes such numerical computa- 
tions affordable also when dynamical quarks are included [23]. We come back to these 
details in Sect. [3l 

The remaining factors in eq. (|2,ip describe the dependence on L. They are called 
step scaling functions. In their original version [24], they depended only on Lj (or 
equivalently a renormalized coupling g(Li)), but here we have an additional dependence 
on the mass of the heavy quark. It is convenient to replace the latter by the dimensionless 
observable 

x = \ \ _|_ q ( i \ (2 2) 

Lmp S (mh,L) Lm^ \( L m h )' 2 J ' 

constructed from a finite volume pseudoscalar heavy-light mass, mps- It will then also 
be used as the HQET expansion parameter instead of the inverse of the heavy quark 
mass. The indicated HQET expansion assumes that L is kept fixed; see e.g. |17|5j for 
more details; it will be used later. First we define the generic step scaling function 

Q(m h ,L) 

a ° {x > L) = 0(m h ,L/s)> (2 - 3) 

(with x from eq. ()2.2j) ) where the scale factor s as well as any other quark masses are 
kept fixed and are not indicated explicitly. 
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In particular, the step scaling function of the pseudoscalar mass itself, 

m PS (m h ,L) 

cr m (x, L) = —r , (2.4 

m PS (m h ,L/s) 

is of central importance. Starting from the experimentally determined mass, m-Q s = 
5.3675(18) GeV and Ljy large enough, it serves to locate the physical points re, via 

x N = l/(L N m Bs ) , Xi_i = s a m (xi, Li) Xi . (2.5) 

The numerical results of all step scaling functions have to be evaluated at these points, 
eq. (|2.ip is then rewritten as 

0(m h ,Loo) = O(m h ,L )a o (x 1 ,Li) ■ ■ ■ a (x N ,L N ) . (2.6) 

Increasing i in eq. (|2.5p successively, the computation of the step scaling functions in the 
relativistic theory at the physical mass requires lattice resolutions Li/ a which become 
larger by a factor s in each step. This is not affordable in practice. Thus the idea of 
|21|22j was to compute ao(x,L{) for a range of x (and thus quark masses) such that 
x > s l Xi ~ xq and to extrapolate x — * x%. In other words in each step (starting from 
Lq) the maximal quark mass which is simulated is reduced by about a factor s. As 
expected from the fact that everywhere one is in the situation x <C 1, the slopes in 
these extrapolations turned out to be rather small and the extrapolations could thus 
be carried out. For an illustration of the x-dependence as it comes out in practice, one 
may look ahead at our final results, Fig. [3l 

Our main point in this paper is that these extrapolations can be turned into inter- 
polations by computing the limiting behavior for small x directly in HQET, 

ao(x,L) =aP(L) + 0(x). (2.7) 

As usual in QCD, this large mass expansion is accompanied by logarithms due to anoma- 
lous dimensions in the effective theory; 0(x) thus stands for at least one power of x 
accompanied by powers of log(x). For cr m , the lowest order term is predicted by the the- 
ory to be one, while for the first order in x, a computation in the static approximation 
of HQET (£ stat = V'h-CoV'h) is required, 

a m (x,L) = l + a s ^ t (L)x + 0(x 2 ). (2.8) 

The static term, which comes from the 0(1/ '(Lm^) 2 ) term in eq. (|2.2p . is not accompa- 
nied by logarithms, see Sect. 13.21 For the case of the decay constant already the lowest 
order term is given by a non-perturbative computation in the static approximation, 

a f (0) (L) = cr f stat (L). (2.9) 

As a further application of this method we compute the mass of the b-quark starting 
from the physical meson mass ?tib s - To this end we define the ratio 

/<«,£) = =£5^ (2.10) 
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Figure 1: The boundary to axial-current correlator /a and the boundary to boundary correlator /i. 
Dirichlet boundary conditions are imposed at xq = and xq = T. 



of the meson mass to the renormalization group invariant (RGI) quark mass, Mh (see 
e.g. [25j for its definition). It provides the connection 



Mb 



p{xo,L )a m (x ll L 1 ) ... a m (x N ,L N ) ' 



(2.11) 



between the physics input t?ib s and the RGI b-quark mass. 

Note that the only approximation made in the above equations is to neglect finite 
size effects in the volume of linear extent L^. 

3 Finite volume observables 



3.1 Relativistic QCD 



Suitable finite volume observables are defined in the QCD Schrddinger functional |26|27j 
with a space-time topology L 3 x T, where T = 2L and C = C = is chosen for the 
boundary gauge fields, and 6 = for the phase in the spatial quark boundary conditions. 

The 0(a)-improved |28I29I3UI31I32] heavy-light correlation functions /a(*)>/p(*) 
and fi are defined and renormalized as in [21j . They are illustrated in Fig. [TJ They 
allow to define a finite volume pseudoscalar meson decay constant and mass via |33|34j 



Fps(m h ,L) 



/a(T/2) 



^ L 3 m PS (m h , L) vTT 



mp S (m h ,L) = — [ln(/ A (T/2 
la 



a))-Hf A (T/2 + a))}. 



(3.1) 
(3.2) 



We remind the reader that we have a fixed ratio T/L = 2. Therefore the time separation 
in the correlation functions grows when L grows. Indeed, as discussed in detail in 
these quantities approach the physical ones in the large L limit, 



lim m PS (m b ,L) = m Bs , 

L— >oo 



lim F PS (m b ,L) = F Bs 

L— >oo 



(3.3) 
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with corrections which (asymptotically) are exponentially small in L. The associated 
step scaling functions are defined as 



, j \ F PS (m h , L) yWsK L) 1 

erf (x, LJ = , x = — , (3.4) 

F PS (m h , L/s) Vm PS (m h , L/s) Lm PS (m h , L) 



and a m as in eq. (|2.4p . 



3.2 HQET 



In the static approximation of HQET, unrenormalized correlation functions /Jj and 
ystat are (j e g ne( j j n complete analogy to the relativistic ones [17] (see [5] for an intro- 
duction). As in these references, we use the RGI static axial current, related to the bare 
one by a factor . It serves to define the RGI ratio , 



W^) = «Gi^=^, (3-5) 

which is related to the QCD decay constant Fps via 

F PS (m h , L)V£ 3 m PS (m h ,L) = -2C PS (M h /A m ) x Y RGl (L) + 0(x) . (3.6) 

The function Cps(Mh/Aj^g), discussed in |35|36j . originates from the matching of QCD 
and the effective theory. In its numerical evaluation we use the anomalous dimension, 
7 PS in the notation of [3B]. With the 3- loop term extracted from [37] its uncertainty 
is estimated to be negligible [36J compared to our other errors. Just like Z^^ Gl , it is 
needed only for Fpg(mh, Lq); it cancels out in the step scaling functions. 
The pseudoscalar finite volume mass has an HQET expansion [17] 

mp S (m h , L) =m h + 5m + T stat (L) + 0(l/m h ) , (3.7) 

with 

T stat (L) = ^[ln(/l tat (T/2 - a)) - ln(/I tat (T/2 + a))] , (3.8) 

where again we do not need to specify the renormalization scheme for m^, but it is 
important that the counterterm dm cancels the linear divergence in r sta t an d the com- 
bination 5m + r sta t(£) is of order Aqcd- Inserting eqs. (|3.6l3.7p into eqs. (|3.4|2.4j) . we 
arrive at eqs. (|2,9l2.8p with the static step scaling functions 

^ {L) = ^W)' (3 - 9) 
aT(L) = L[r stat (L) -r stat (L/s)]. (3.10) 

Here the renormalizations Z^^ Gl x Cps and 5m cancel, which shows that these static 
step scaling functions are not accompanied by any logarithmic terms (in or x). 

In our numerical investigation we will compute them precisely by using the static 
action denoted by HYP2 in [38] (see also [39]), and the corresponding 0(a)-improvement 
coefficients for the static axial current. 
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4 Results in the quenched approximation 



We employ the non-perturbatively 0(a)-improved Wilson action |29|31j . The data at 
finite heavy quark mass are taken from |21|22j . As there, we choose N = 2 steps, s = 2 
and Lq = 0.4 fm. The length scale is set by ro = 0.5 fm [30] using the parametrizations 
of ro I a as a function of the bare coupling go from |41|42j . The light quark mass is set to 
the strange quark mass by fixing the RGI-mass to M s = 0.1346(55) GeV as previously 
determined from the Kaon mass in the quenched approximation [43]. The RGI-mass 
is related to the bare one by a non-perturbatively computed renormalization factor 
Z m [25], see e.g. [21] for details. 



L[fm] 


X 


<r m (x,L) 


(Tf(x,L) 


1.6 


0.0581 
0.0670 
0.0720 


1.069(5) 
1.081(6) 
1.087(7) 


0.929(32) 
0.912(27) 
0.900(24) 


0.8 


0.0804 
0.0884 
0.1204 


1.012(6) 
1.014(6) 
1.018(8) 


0.4198(45) 
0.4193(45) 
0.4169(43) 






p(x,L) 


<p(x, L) 


0.4 


0.0933 
0.0990 
0.1472 
0.2768 
0.2885 


0.744(09) 
0.754(09) 
0.837(12) 


3.120(45) 
3.097(45) 
2.911(43) 
2.534(40) 
2.505(40) 



Table 1: Finite mass observables after continuum extrapolation. Physical units are set through 
ro = 0.5 fm. Statistical errors of x due to tops have been changed to errors in the x-dependcnt 
observables. 

4.1 At finite heavy quark mass 

The data of |21|22j have been reanalyzed. The step scaling functions were first defined 
at a fixed value of ro-Mh as in those references. Their continuum limit was taken by 
an extrapolation linear in (a/L) 2 , making use of different definitions of at finite 
lattice spacing and of the fact that the continuum limit is independent of such details. 
Correlations between observables computed on the same gauge configurations were taken 
into account. The statistical uncertainties of the regularization dependent part of the 
renormalization constants and the lattice spacing were included before performing the 
continuum limit extrapolations, the uncertainty of the regularization independent part 
of the renormalization constants is added in the continuum limit; all these do not appear 
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as a separate uncertainties, rather they are included in the quoted errors. For their 
detailed accumulation we refer to [H]. An impression on the quality of the continuum 
extrapolations is easily obtained from the graphs in [21|22|44] . Since here our emphasis 
is on the use of the static approximation, we do not reproduce those details. The 
continuum values of the step scaling functions were then interpolated in the pseudoscalar 
mass to a few selected values of x. These are listed in Table [1] together with p, eq. (|2.10p 
and 

<p{x,L) =L 3 / 2 F PS y^. (4.1) 



4.2 In static approximation 

We turn to the main new element in our numerical computations. We start with 
the static step scaling function <7^ at (L2)) requiring the computation of r sta t(^i) and 
rstat(-^2)) for several fixed values of go followed by a continuum extrapolation. How- 
ever, it is a central element of our strategy that L 2 ~ 1.6 fm is large enough such that 
finite volume effects are negligible. Thus we can replace r sta t(-^2) by E stll t, the "mass" of 
a static-strange bound state in large volume which is known from |19|45|46] in the range 
6.0219 < P = Q/gl < 6.4956. We have computed r stat (Li) for L/a = 8,10,12,16,24, 
spanning a wider range in (3 and allowing easily for an interpolation to the values of 
(3 where S s tat is known. All of this was done for the HYP2 static action [38] and for 
the tree- level as well as the one- loop improved static-light axial current. Differences 
between the two turned out to be far below our statistical precision of order 1 — 3MeV. 
The continuum extrapolation a^ at (L 2 ) , listed in Table [2j is well controlled, see Fig. [2j 
Similarly we profit from previous work in large volume in the computation of 
<T£ tat (L2), eq. (13.91) . In that case the continuum value 

Yrgi{L 2 ) = -4.65(19) , (4.2) 

is known from [45|46j . It remains to compute 

W^i) = -^j L TxZ^ t (L'/a,g )xY(L 1 /a,g ), L' = l/fi. (4.3) 

Here 

nL/a , m] = .^n ( , 4) 

is the unrenormalized version of eq. (|3.5p and Z^ tat (L'/a, go) is the factor, introduced in 
|35j . to renormalize the static axial current in the ("new") SF scheme, non-perturbatively 



at renormalization scale fi = 1/L' . Finally re ^ es an y matrix element of the ax- 

ial current in the SF scheme to the RGI matrix element. At the renormalization point 
L' = L\ its non-perturbative value 
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Li/a 


a 


<zr stat (Li) 


ar stat (i 2 ) 


< & \L 2 ) 


8 


5.9598 


0.3183(8) 






8.92 


6.0219 


0.3000(5) 


0.4053(49) 


1.878(88) 


10 


6.0914 


0.2805(6) 






12 


6.2110 


0.2533(6) 






13.41 


6.2885 


0.2360(5) 


0.3011(33) 


1.745(89) 


16 


6.4200 


0.2114(7) 






16.64 


6.4500 


0.2067(6) 


0.2564(09) 


1.654(35) 


17.64 


6.4956 


0.1997(6) 


0.2461(14) 


1.637(54) 


24 


6.7370 


0.1722(16) 






continuum 






1.561(53) 



Table 2: Numerical results in static approximation for L\ = 0.8 fm, L 2 = 2L\. The rows with 
non-integer L\ja list interpolated values for a,r stat (Li), while ar stat (L2) are the large volume 
numbers of [19l45|46j . 



Lo/a 13 


a r stat (L ) 


ar stat (Li) 


< at (Li) 


Y(L ) 


Y{L t ) 


CT f stat (il) 


6 6.2110 


0.2272(9) 


0.2558(18) 


0.343(24) 


-1.805(03) 


-2.221(13) 


0.4350(27) 


8 6.4200 


0.1958(9) 


0.2154(11) 


0.315(23) 


-1.837(05) 


-2.266(13) 


0.4361(26) 


12 6.7370 


0.1561(8) 


0.1663(17) 


0.245(46) 


-1.881(06) 


-2.279(28) 


0.4284(55) 


16 6.9630 


0.1355(7) 


0.1426(14) 


0.230(50) 


-1.899(07) 


-2.344(35) 


0.4366(67) 


24 7.3000 








-1.918(10) 






continuum 






0.233(36) 






0.4337(44) 



Table 3: Numerical results in static approximation for Lq = 0.4 fm and L\ = 2Lq. 



La 
a 





ysta.tr L a \ 


O 


& 


^I tat (^,3o) 


8 


6.4200 


0.8745(21) 


12 


6.2110 


0.7904(38) 


12 


6.7370 


0.8534(10) 


16 


6.4200 


0.7672(45) 


16 


6.9630 


0.8408(21) 


24 


6.7370 


0.7651(53) 


24 


7.3000 


0.8308(21) 


32 


6.9630 


0.7556(48) 



Table 4: Renormalization factors for the static axial current at renormalization scales fi = I/Lq 
and n — l/Li with Lq = 0.4 fm and L\ = 0.8 fm. 
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0.01 0.02 0.01 0.02 

(a/L)» (a/L)» 

Figure 2: Continuum extrapolations of static results. The extrapolated values with their errors are 
shown at a/L = 0. 



is easily extracted from the results in [3S]. We have computed the missing factors 
Zj£ at (L'/a,go),Y(Li/a,go) for various values of Li/a, setting L' = L\, see Table |U 
Note that following the exact definition of [35] , = 1/2 , T = V is employed for Z*F 
and the computation is carried out at zero (light) quark mass - in contrast to the 
evaluation of Y (and all other quantities). 
The continuum limit (Fig. [2]) 

Yrgi(Li) = -1.628(19) (4.6) 

is combined with (|4,2p to get 

a f stat (L 2 ) = 1.010(43). (4.7) 

In the computation of the static step scaling functions (j stat (Li) (Table [3l Fig. [2]) 
we followed straightforwardly their definitions. Finally, 

<I>K<;) 0.846(6). (4.8) 



*sp(1/£o) 

from [55] together with Y(L /a, g ) , Z^ at (L /a ; go), Tables [3] and [H yields 

Yr G i(L ) = -1.347(13) (4.9) 
by a continuum extrapolation again illustrated in Fig. [2l 
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4.3 Interpolation to the physical point 

We now combine the static results with the relativistic ones, through linear and quadratic 
interpolations in x. Namely we fit for the parameters rrij(Li) and ej(Li) in 

1 + mi (Li) x + m 2 (Li) x 2 , (4.10) 
mi (L f ), (4.11) 
e (L) + ei(Li) x + e 2 {U) x 2 , (4.12) 

and then insert the fit functions eq. (|4.10p and eq. (|4.12p into eq. (|2.1ip and eq. (|4.13p . 
Note that the first two equations are fit together. The static <7^ at (Lj) enter eq. (|4.1ip as 
data points and af' a,t {Li) are data at x = in eq. (|4.12p . As seen in Fig. [21 the quadratic 
terms are moderate in the whole range and in particular at the physical points xi the 
differences between the static results and the interpolated ones are rather small. As an 
illustration of the effect of the static results we also carry out an analysis where they 
are not taken into account. The numbers in Table [5] show that the statistical errors in 
the step scaling functions are significantly reduced by including the static constraints. 
Furthermore we can perform the consistency check of including quadratic terms only 
when the static constraints are used. The agreement between linear and quadratic 
interpolations is very reassuring. 



i 


cr m (x 


it Li) 


o-{{xi, 


Li) 


Fit 


2 


1.0330(11) 


1.0258(21) 


0.985(31) 




quadratic 


2 


1.0319(11) 


1.0276(22) 


0.911(29) 


1.002(54) 


linear 


1 


1.0092(18) 


1.0074(33) 


0.4243(36) 




quadratic 


1 


1.0093(15) 


1.0072(32) 


04260(31) 


0.4223(48) 


linear 



Table 5: Step scaling functions inter/extra-polated to the physical points x-i — 0.022974(8), 
xi = 0.04746(5) and x^ CD = 0.04741(10), where the latter originates from the fits to only the 
finite heavy quark mass data. The left-side column of each a is the number including the static 
constraint, the right-side one is without. 

For p(x,Lq) and Lpg(x,Lo) the relativistic simulations straddle the physical point 
x = xq and, for the decay constant, the static data do not sensitively improve the 
precision on the interpolated point. However, as an illustration how HQET does describe 
these quantities, we also show eq. (|4.9p together with the data at finite x in Fig. [3J in 
that case the interpolation yields L Q Lpg- v /mpg/(2Cps) = 1.279(17) or c/?(xo,Lo) = 
3.107(41) and p{x ,L ) = 0.7485(9). 

Our final large volume results from eq. (|2.1ip and 

L Bs = ip(x , L ) cr f (xi, Li) <7f (x 2 , L 2 ) L 3/2 m^J 2 (4.13) 



0~m(x,Li) — 

<\U) = 
a f (x,Li) = 
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Figure 3: Interpolations to the physical points are shown by the filled circles. For cr m , the static 
constraints are illustrated as the error band of the static result 1 + fr^ at i. On the right hand side, the 
static results enter as data points at x — 0. 



are 

F Bs = 191(6) MeV , M b = 6.88(10) GeV => m b (m b ) = 4.42(6) GeV . (4.14) 

Here the conversion to the running mass in the MS-scheme is done with the 4-loop RG 
equations (for N f = and A|| = 238(19) MeV [25]). 

5 Conclusions and outlook 

We have followed a general strategy for computing B-meson observables. Starting from 
a finite volume, where the observables are straightforwardly computable in relativistic 
lattice QCD, we evaluated step scaling functions which describe the finite size effects. 
The latter are not directly computable at the physical points since for accessible lattices 
am b > 1. Previously these functions have either been computed by an extrapolation in 
the heavy quark mass to the physical m b |21|22j or they have been computed in HQET 
|17|19j . Here we have demonstrated how the two approaches can be combined to further 
increase precision and confidence in the results. 

Figure El which is a continuum graph, shows that the static (lowest order HQET) 
results match very well onto the finite mass step scaling functions. We therefore have 
excellent control over the heavy quark mass dependence - if desired from below the 
charm quark mass to the b-quark mass and beyond. 
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Our final numbers for decay constant and b-quark mass, eq. (|4.14p . agree well with 
the previous estimates of |21|22I1 9 45 46] where the same experimental data was used 
as input H 

In our results, Fig. El one notices that the corrections to the static approximation 
are very small at the b-quark mass. This represents an intriguing demonstration of 
the precision and usefulness of HQET for B-physics. Although our exercise was in the 
quenched approximation, such a qualitative result may well be carried over to (full) 
QCD. 

Concerning the application of the strategy to QCD, the attentive reader will have 
noticed that in our computations we extensively relied on the knowledge of a reference 
scale {ro/a) over a large range of lattice spacings a. This luxury is not available in 
full QCD - and will not be for a while to come. However, with the knowledge of the 
running coupling of [47j . one can properly set the scale also for small lattice spacings. 
We further note that the finite volume computations which are needed in this strategy 
require a significantly smaller effort than the large volume ones. 

We therefore conclude that the here investigated method is very promising for the 
near future where we expect that high precision can be reached for B-physics. Note 
that the strategy may be extended to other observables such as mass splittings |48|49j 
and form factors [50J. 
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